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Abstract 

We make an observation about Galilean transformation on a 1-D mass 
H ' variable systems which leads us to the right way to deal with these systems. 

■ - - Then using this observation, we study two-bodies gravitational problem where 

the mass of one of the bodies varies and suffers a damping-anti damping effect 
due to star wind during its motion, for this system, a constant of motion, 
a Lagrangian and a Hamiltonian are given for the radial motion, and the 
period of the body is studied using the constant of motion of the system. 
An application to the comet motion is given, using the comet Halley as an 
example. 
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1 Introduction 

There is not doubt that mass variable systems have been relevant since the founda- 
tion of the classical mechanics and modern physics (Lopez et al 2004). These type of 
systems have been known as Gylden-Meshcherskii problems (Gylden 1884; Meshch- 
erskii 1893, 1902; Lovett 1902; Jeans 1924; Berkovich 1981; Bekov 1989; Prieto and 
Docobo 1997), and among these type of systems one could mention: the motion 
of rockets (Sommerfeld 1964), the kinetic theory of dusty plasma (Zagorodny et al 
2000), propagation of electromagnetic waves in a dispersive-nonlinear media (Seri- 
maa et al 1986), neutrinos mass oscillations (Bethe 1986; Commins and Bucksbaum 
1983), black holes formation (Helhl et al 1998), and comets interacting with solar 
wind (Daly 1989). This last system belong to the so called "gravitational two-bodies 
problem" which is one of the most studied and well known system in classical me- 
chanics (Goldstein 1950). In this type of system, one assumes normally that the 
masses of the bodies are fixed and unchanged during the dynamical motion. How- 
ever, when one is dealing with comets, beside to consider its mass variation due to 
the interaction with the solar wind, one would like to have an estimation of the the 
effect of the solar wind pressure on the comet motion. This pressure may produces a 
dissipative-antidissipative effect on its motion. The dissipation effect must be felt by 
the comet when this one is approaching to the sun (or star), and the antidissipation 
effect must be felt by the comet when this one is moving away from the sun. To deal 
with these type of mass variation problem, it has been proposed that the Newton 
equation must be modified (Sommerfeld 1964; Plastino and Muzzio 1992) since the 
system becomes non invariant under change of inertial systems (Galileo transforma- 
tion). In this paper, we will make first an observation about this statement which 
indicates the such a proposed modification of Newton's equation has some problems 
and rather the use of the original Newton equation is the right approach to deal 
with mass variation systems, which it was used in previous paper (Lopez 2007) to 
study two-bodies gravitational problem with mass variation in one of them, where 
we were interested in the difference of the trajectories in the spaces {x, v) and {x,p). 
As a consequence, there is an indication that mass variation problems must be dealt 
as non invariant under Galilean transformation. 



Secondly, , we study the two-bodies gravitational problem taking into consider- 
ation the mass variation of one of them and its damping-anti damping effect due 
to the solar wind. The mass of the other body is assumed big and fixed , and the 
reference system of motion is chosen just in this body. In addition, we will assume 



that the mass lost is expelled from the body radially to its motion. Doing this, 
the three-dimensional two-bodies problem is reduced to a one-dimensional problem. 
Then, a constant of motion, the Lagrangian, and the Hamiltonian are deduced for 
this one-dimensional problem, where a radial dissipative-antidissipative force pro- 
portional to the velocity square is chosen. A model for the mass variation is given, 
and the damping-anti damping effect is studied on the period of the trajectories, 
the trajectories themselves, and the aphelion distance of a comet. We use the pa- 
rameters associated to comet Halley to illustrate the application of our results. 



2 Mass variation problem and Galileo transfor- 
mation. 

To simplify our discussion and without losing generality, we will restrict myself to 
one degree of freedom,. Newton equation of motion is given by 

-{m{t)v) = F{x,v,t), (1) 

where m(t) v is the quantity of movement, F is the total external force acting on the 
object, m{t) and v = dx/dt are its time depending mass and velocity of the body 
(motion of the mass lost is not considered). Galileo transformations to another 
inertial frame (S") which is moving with a constant velocity u respect our original 
frame S are defined as 

X = X — ut (2a) 

t' = t (2b) 

which implies the following relation between the velocity seen in the reference system 
S, V, and the velocity seen in the reference system 5", v', 

v' = V — u. (3) 

Multiplying the last term by m{t') and making the differentiation with respect to 
t', one gets 

A(m(t>')=FV,<0, (4) 



where F' is given by 

F\x\ v', t') = F{x' + ut', v' + u, t') - «^^. (5) 

Therefore, Eq. 1 and Eq. 4 have the same form but the force is different since 
in addition to the transformed force term F{x' + ut',v' + u,t'), one has the term 
udmit') / dt' . This non invariant form of the force under Gahlean transformation has 
lead to propose (Sommerfeld 1964; Plastino and Muzzio 1992) that Newton equation 
(1) to modify Newton's equation of motion for mass variation objects, to keep the 
principle of invariance of equation under Galilean transformations, of the form 

m{t)'^ = F{x,v,t) + w^^ , (6) 

where w is the relative velocity of the escaping mass with respect the center of mass 
of the object. When one does a Galilean transformation on this equation, one gets 

m{t')^^ = F\x\v\t'), (7) 

where F' is given by 

F'(x', v\ t') = F{x' + ut\ v' + n, t') + w;^^, (8) 

which has the same form as Eq. 6. However, assume for the moment that w = 
constant and F = 0. So, from Eq. 6, it follows that 

v{t)=v, + \n(^^]\ (9) 

where niQ = ?7i(0). In this way, if we have a mass variation of the for ni(t) = ttiqC""* 
(for example), one would have a velocity behavior like 

v(t) = vq — wat, (10) 

which is not acceptable since one can have w > 0, w = and v < depending on 
the value wat. Even more, since for F = 0, the equation resulting in the reference 
system S' is the same, i.e. in S' one gets the same type of solution, 

v'{t') = Vo + \n('^^) (11) 



which is independent on the relative motion of the reference frames, and this must 
not be possible due to relation (3). 

In addition, it worths to mention that special theory of relativity can be seen as 
the motion of mass variation problem, where the mass depends on the velocity of 
the particle of the form m{v) = mo(l — t'^/c2)~^/^, with c being the speed of light. 
This system is obviously not invariant under Galilean transformation, and given 
the force, Newton's equation motion is always kept in the same form to solve a 
relativistic problem, d{m{v)v) / dt = F{x,v,t), (C M0ller 1952, Lopez et al 2004). 

3 Mass variation and equations of motion. 

Having explained and clarify the problem of mass variation, Newton's equations 
of motion for two bodies interacting gravitationally, seen from arbitrary inertial 
reference system, and with radial dissipative-antidissipative force acting in one of 
them are given by 

d ( dYx\ Gmim2 

'dt 



Jt^"^'- 



and 
d 



dT2 

dt V^'^~dt 



Gmim2 

|r2 -ri|3 



r2 - ri 1 



ri -r2 



7 
ki -r2| 



r2j 



(12a) 



d|ri — r2| 
dt 



r2 - ri 



(12b) 



where rrii and 7722 are the masses of the two bodies, ri = (xi,?/i, Zi) and 
r2 = {x2,y2,Z2) are their vectors positions from the reference system, G is the 
gravitational constant {G = 6.67 x 10~^^m^/Kg s^), 7 is the nonnegative constant 
parameter of the dissipative-antidissipative force, and 



^2 



^2 



Vi^2 - Xiy + (?/2 - yiy + {Z2 - ZiY 



is the Euclidean distance between the two bodies. Note that if 7 > and 
d\Y\ — r2|/dt > one has dissipation since the force acts against the motion of the 
body, and for d\T\ — r2|/dt < one has anti-dissipation since the force pushes the 
body. If 7 < this scheme is reversed and corresponds to our actual situation with 
the comet mass lost. 



It will be assumed the mass mi of the first body is constant and that the mass 
m2 of the second body varies. Now, It is clear that the usual relative, r, and center 



of mass, R, coordinates defined as r = r2 — ri and R = {niir i+ 7712^2) / {fni + 1712) are 
not so good to describe tlie dynamics of this system. However, one can consider the 
case for mi ^ m2 (which is the case star-comet), and consider to put our reference 
system just on the first body (ri = 0). In this case, Eq. (12a) and Eq. (12b) are 
reduced to the equation 



"^2377 = ^ — r - msr - 7 



dr'^ 



dt 



(13) 



where one has made the definition r = r2 = {x,y,z), r is its magnitude, r = 
a/x^ + y"^ + z"^, and r = r/r is the unitary radial vector. Using spherical coordinates 

x = rsm6cos(f, y = rsinOsirup , z = rcos6, (14) 

one obtains the following coupled equations 



/■■ n2 -2 ■ 2n\ Gmim2 . . .2 /i c; ^ 

m2[r — rU — np sm u) = ^ m2r — ^yr , (15a) 



m2{2r9 + rO — np^ sin 6 cos 6) = —7fi2r9 , (15b) 

and 

m2{2ripsin6 + np sin 6 + 2rip6 cos6) = —7h2ripsin6 . (15c) 

Taking yj = as solution of this last equation, the resulting equations are 

... A2N Gmim2 . . .2 ,^„ s 

fji2{r — ru ) = m2r — ^yr , (loaj 

and 

m2{2re + rO) + m2r6 = . (16b) 

From this last expression, one gets the following constant of motion (usual angular 
momentum of the system) 

k = m2r^e , (17) 

and with this constant of motion substituted in Eq. (16a), one obtains the following 
one-dimensional equation of motion for the radial part 

d'^r Grrii m2 f dr\ 7 .2 , ^0 n on 



dt'^ r^ 777,2 \dt J r?72 m^r^ 



Now, let us assume that ?7i2 is a function of the distance between the first and the 
second body, m2 = rra2(r). Therefore, it follows that 

7712 = "Tijr , (19) 

where m'g is defined as rra'g = dm2/dr. Thus, Eq. (18) is written as 
dV Gmi , Ij m'2 + 7 /rfr \ ^ 



(it^ r^ ?72|r'^ r7i2 yc/t 

which, in turns, can be written as the following autonomous djTiamical system 

dr dv Grrii ll m, + 7 , 

w7=^' w7 = - + -h-— -"" ■ (21) 

dt dt r'^ TTigr-* m2 

Note from this equation that m^ is always a non-positive function of r since it 
represents the mass lost rate. On the other hand, 7 is a negative parameter in our 
case. 

4 Constant of Motion, Lagrangian and Hamilto- 
nian 

A constant of motion for the dynamical system (21) is a function K = K{r, v) which 
satisfies the partial differential equation (Lopez 1999) 



OK 

dr 



-Gmi ll _ m'^ + 'j 2 



'" 0. (22) 



dv 

The general solution of this equation is given by (John 1974) 

K{x,v) = F{c{r,v)), (23) 

where F is an arbitrary function of the characteristic curve c{r, v) which has the 
following expression 

c{r, v) = m2(r)e2^"W^;2 ^ j f^Gm^ _ _^_\ ^l^r)e^iKr)dr^ (24) 

and the function A(r) has been defined as 

Kr)=f^ (25) 

J m2(r) 

7 



During a cycle of oscillation, the function m2{r) can be different for the comet 
approaching the sun and for the comet moving away from the sun. Let us denote 
m2+{'^) for the first case and ?7i2_(r) for the second case. Therefore, one has two 
cases to consider in Eq. (23) which will denoted by (±). Now, if m2± denotes the 
mass at aphelium (+) or perielium (-) of the comet, F{c) = c^/2m2^ represents the 
functionality in Eq. (23) such that for 7712 constant and 7 equal zero, this constant 
of motion is the usual gravitational energy. Thus, the constant of motion can be 
chosen as K^ = c(r,t') /2?7i2±, that is. 



K^ = '^e^''^'''-' + Klfir\ (26a) 

where the effective potential Veff has been defined as 

This effective potential has an extreme at the point r^, defined by the relation 

nmlin) = -^ (27) 

Grrii 

which is independent on the parameter 7 and depends on the behavior of 7722 (r). 
This extreme point is a minimum of the effective potential since one has 

-Jl] > 0. (28) 

Using the known expression (Kobussen 1979; Leubner 1981; Lopez 1996) for the 
Lagrangian in terms of the constant of motion, 

^ / N f Kir, v) dv , , 

L{r,v)=vj "-y , (29) 

the Lagrangian, generalized linear momentum and the Hamiltonian are given by 



''2± 



^Ljr) V ^2,X+(r) 



p='^:^^^L^e^^>^±ir)^ (31) 
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and 



,o ^2 



^^-aTH-"""^"-""- 



(32) 



•'2±\ 



The trajectories in the space (x, v) are determined by the constant of motion (26a). 
Given the initial condition {to^Vo)-, the constant of motion has the specific value 



and the trajectory in the space (r, v) is given by 



(33) 



f = ±< 



"^i±(^) 



-7A±(r-) 



^o - Ktfir) 



1/2 



(34) 



Note that one needs to specify 9o also to determine Eq. (17). In addition, one 
normally wants to know the trajectory in the real space, that is, the acknowledgment 
of r = r{6). Since one has that v = dr/dt = {dr/d6)6 and Eqs. (17) and (34), it 
follows that 



eir) = 9o + 






m2± 



(^)e7A±{r)^^ 



\/2mf± Jr„ t^Jk^ - Vtff{r\ 



The half-time period (going from aphelion to perihelion 
be deduced from Eq. (34) as 



^1/2 



v/2^ 



2± "'r-i 



m2± 



V)eT^±('')rfr 



K^ - KU^) 



(35) 
or backward (-)) can 

(36) 



where ri and r2 are the two return points resulting from the solution of the following 
equation 

V^%{r,) = K^ ^ = 1,2. (37) 

On the other hand, the trajectory in the space {r,p) is determine by the Hamiltonian 
(32), and given the same initial conditions, the initial po and H^ are obtained from 
Eqs. (32) and (31). Thus, this trajectory is given by 



p = ±w'^^^^^e-^-(^) 



m 



2± 



H^ - Klfir) 



1/2 



(3^ 



It is clear just by looking the expressions (34) and (38) that the trajectories in the 
spaces (r, t>) and {r,p) must be different due to comphcated relation (31) between v 
and p (Lopez 2007). 



5 Mass- Variable Model and Results 

As a possible application, consider that a comet looses material as a result of the 
interaction with star wind in the following way (for one cycle of oscillation) 



m2±(r) 



m2-(r2(i_i), 



'm2+[r2i^i, 



l-e- 



incom,ing{+) v < 



6(1-6 "^'' "'zi-i) j outgoing{-) f > 



(39) 



where the parameters 6 > and a > can be chosen to math the mass loss rate in 
the incoming and outgoing cases. The index "i" represent the it h- semi- cycle, being 
^2(j-i) and r2i_i the aphelion(ra) and perihelion(rp) points [tq is given by the initial 
conditions, and one has that m2-{ro) = rrio). For this case, the functions A+(r) and 
A_(r) are given by 

A+(r) = -^ In ( e"" - 1 ) , (40a) 



arUa 



and 



A_(r) 



-1 



ar 



+ In (rrip — 6(1 — e 



-a{r-rp) 



)) 



(40b) 



a{b — Trip) 

where we have defined rria = m2{ra) and rrip = m2{rp). Using the Taylor expansion, 
one gets 



and 



g27A+(r) _ ^2-fr/ma 



,2'yX-ir) 



arUa 



2"fr 



g (b-mp) 



2 aTTia yarria 

27 



(41a) 



27 



{nip — 6)'»('"P-'') 
1 27 



-a{r-rp) 



a{mp — b) nip 



27 



2 a{nip — b) \a{mp — b) 
10 



1 



g-2a(r-r-p) 

{nip - by 



(41b) 



The effective potential for the incoming comet can be written as 



K'ffir) 



Gmirria Ig 1 



r 2ma r^ 

and for the outgoing comet as 

Grriima tl 1 



e^^'l"^- + Wi{-i,a,r), (42a) 



K7/(0 



2mn r^ 



27r 
(mp—b) 



(D{mp 



27 



{mp - fe)"!™?- 



+ iy2(7,«,0, (42b) 



b) 



where Wi and 14^2 are given in the appendix A. 

We will use the data corresponding to the sun mass (1.9891 x lO^'^Kg) and the 
Halley comet ( Cevolani et al 1987, Brandy 1982, Jewitt 2002) 

rric ^ 2.3 x lO^^Kg, Vp ^ 0.6 au, r^ ^ 35 au, k ^ 10.83 x lO"^^ Kg ■ m'^/s, (34) 

with a mass lost of about 6m ^ 2.8 x lO^^Kg per cycle of oscillation. Although, 
the behavior of Halley comet seem to be chaotic (Chirikov and Vecheslavov 1989), 
but we will neglect this fine detail here. Now, the parameters a and "6" appearing 
on the mass lost model, Eq. (39), are determined by the chosen mass lost of the 
comet during the approaching to the sun and during the moving away from the sun 
(we have assumed the same mass lost in each half of the cycle of oscillation of the 
comet around the sun). Using Eq. (42a) and Eq. (42b) in the expression (34), the 
trajectories can be calculated in the spaces (r, v) . Fig. 1 shows these trajectories 
using 6m = 2 x 10^^ Kg (or 6m/ m = 0.0087%) for 7 = and (continuos line), and 
for 7 = —3 Kg/m, (dashed line), starting both cases from the same aphelion dis- 
tance. As one can see on the minimum, dissipation causes to reduce a little bit the 
velocity of the comet , and the antidissipation increases the comet velocity, reaching 
a further away aphelion point. Also, when only mass lost is considered (7 = 0) the 
comet returns to aphelion point a little further away from the initial one during 
the cycle of oscillation. Something related with this effect is the change of period 
as a function of mass lost (7 = 0). This can be see on Fig. 2, where the period is 
calculated starting always from the same aphelion point (r^). Note that with a mass 
lost of the order 2.8 x lO^^Kg (Halley comet), which correspond to 6m/m, = .12%, 
the comet is well within 75 years period. The variation of the ratio of the change 
of aphelion distance as a function of mass lost (7 = 0) is shown on Fig. 3. On Fig. 
4, the mass lost rate is kept fixed to 6m/m, = 0.0087%, and the variation of the 
period of the comet is calculated as a function of the dissipative-antidissipative pa- 
rameter 7 < (using I7I for convenience). As one can see, antidissipation always 

11 



wins to dissipation, bringing about the increasing of the period as a function of 
this parameter. The reason seems to be that the antidissipation acts on the comet 
when this ones is hghter than when dissipation was acting (dissipation acts when 
the comet approaches to the sun, meanwhile antidissipation acts when the comet 
goes away from the sun). Since the period of Halley comets has not changed much 
during many turns, we can assume that the parameter 7 must vary in the interval 
(—0.01, 0]Kg/m. Finally, Fig. 5 shows the variation, during a cycle of oscillation, of 
the ratio of the new aphelion (r^) to old aphelion (r^) as a function of the parameter 
7- 



6 Conclusions and comments 

We have shown that the proposed modified Newton equation for mass variation 
systems has some problems. Therefore, we have considered that it is better to keep 
Newton's equations of motion for mass variable systems to have a consistent ap- 
proach to these problems. Having this in mind, the Lagrangian, Hamiltonian and 
a constant of motion of the gravitational attraction of two bodies were given when 
one of the bodies has variable mass and the dissipative-antidissipative effect of the 
solar wind is considered. By choosing the reference system in the massive body, 
the system of equations is reduce to 1-D problem. Then, the constant of motion, 
Lagrangian and Hamiltonian were obtained consistently. A model for comet-mass- 
variation was given, and with this model, a study was made of the variation of 
the period of one cycle of oscillation of the comet when there are mass variation 
and dissipation-antidissipation. When mass variation is only considered, the comet 
trajectory is moving away from the sun, the mass lost is reduced as the comet is 
farther away (according to our model), and the period of oscillations becomes bigger. 
When dissipation-antidissipation is added, this former effect becomes higher as the 
parameter 7 becomes higher. 

It is important to mention that if instead of loosing mass the body would had win- 
ning mass, the period of oscillation of the system would decrease. One can imagine, 
for example, a binary stars system where one of the star is winning mass from the 
interstellar space or from the other star companion, and due to this winning mass, 
the period of the star would decrease depending on how much mass the star is ab- 
sorbing. 



12 



13 



7 Appendix A 

Expression for Wi and W2'- 



Wi = — ^i -^-^ '- + apEi{apr) - 2ap{p - l)Ei{{-4 + p)ar) 



m. 



2+ 



2r 



'^2^^^E,{{-A + p)ar) + ^^^-^ [e^^'^)-^ + 3a(l -p)rE,((p - 3)ar)] 



p(p + 3) 



2 

p + 2 



,{p-2)ar 



+ a{p~2)Ei{{p-2)ar) 



^^{p-2)ar ^ ^^p _ iyE^{{p - l)ar)] 



h I P{P ~ 1) Jp-2)ar 



2m2+r2 1 2 



(^-2)"" -pe^P-^)"" - ap(p - l)e(P-2)"^- + '^^ — l^e^"" 



,^M^J)!:e(p-2)- +p«re(P-^)"^^ -p2«re(^-i)"^^ -p2«V2E,(par) 



a^{p — 2Yp{p — l)r' 



Eiijyp — 2)ar) + pa^r'^Ei{{p — l)a;r) 



-2a'^p^r^Ei{{p - l)ar) + p^a^r^Ei[{p - l)ar) 



{Al) 
where ma is the mass of the body at the aphehon, and we have made the definitions 

27 



p 



arrin 



{A2) 



and the function E^ is the exponential integral 

EAz) = 



OOg-f 



dt 

t 



(A3) 
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Wo 



Gml_ \ e^'^-^^' 



TTln 



+qaEi (gar) 



2{mp + ag) m^ + ag 

1^ [.<-." -a.(,-3).B.((,-3)a.)] 



ge 



qar 



{rup + ag)r 



,((j-3)ar 



a{q - 3)rEi{{q - 3)ar)] - 2aEi{{q - 2)ar) 



+ag£;,((g - 2)ar) - ^ -^——E,{{q - 2)ar) 

^ [nip + aqY ^ ' 

^^ ^^^E.((g - 2)ar) E,{{q - 2)ar) 



2(mp + ag)2 
2g2ae«'^'' 



rrtp + ag 



(mp + aq)r 



qe 



qar 



[mp + aq)r 

h 

2m2_^{mp + aqY 



Ei{{q - 2)ar) + - [e^'^"^)"" - (g - l)arEi{{q - l)ar)] 
^^{q-i)ar _ ^^ _ i)arE,{{q - l)ar)] 



qae 



qar 



+ q^a^Ei [qar) 



+ ^:^ ' .. . [-1 + 2ar - g«r + (2 - g)2«V2e(2-'')"'^E.((g - 2)ar)] 






2{mp + aqYr 



{nip + ag)r2 



[-1 + ar + qar + (g - l)2aV2e(i-'?)°'"Ei((g - l)ar)] 



(A4) 
where mp is the mass of the body at the perihehon, and we have made the definition 
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Q 



a{mp — h) 



{Ah) 
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Figure 1: Trajectories in the {r,v) space with 6m/m = 0.009. 
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Figure 2: Period of the comet as a function of the mass lost ratio. 




Figure 3: Ratio of aphelion distance change as a function of the mass lost rate. 




Figure 4: Period of the comet as a function of the parameter 7. 
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Figure 5: Ratio of the aphelion increasing as a function of the parameter 7. 
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